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We have studied tlieoretically tiie Bose-Einstein condensation (BEC) of two-dimensional excitons 
in a ring with small random width variation. We derive a nonlinear Gross-Pitaevkii equation (GPE) 
for such a condensate. Our numerical solution of the ground state of the GPE displays a necklace- 
like structure in the presence of small random variation of ring width. This is a consequence of the 
interplay between random potential and the nonlinear repulsive term of the GPE. Our result suggests 
that the formation of ring and necklace-like structures observed recently in the photoluminescence 
of quantum wells may be stable even in the BEC phase. 
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Recent advances in the Bose-Einstein condensation 
(BEC) 0,01 01 of atoms were made possible by exploring 
the confinement of atomic gases within potential traps 
0. In semiconductor physics, the possibihty of exciton 

I BEC has been explored for a long time 0; 0- the 
dilute limit nag << 1 (a^ is the exciton Bohr radius, D 

[ the dimensionality, and n is the exciton density) excitons 
are expected to undergo the BEC ,5j . Recently, the pos- 
sibility of BEC of excitons in potential traps of semicon- 
ductor quantum wells has been explored [E Hi IS ^3 ■ 
The quantum degeneracy temperature is expected to be 
about IK, much higher than that of atoms due to the 
small mass of the exciton. Butov et al. Q i Snoke et 

] al. and Rapaport et al. ^ have observed a circu- 
lar ring in exciton photoluminescence (PL) from these 
quantum wells, where the radius of the ring (typically 
10 ~ 100/im) depends on the laser excitation intensity. 
A mechanism for the formation of the ring is proposed 
that involves carrier imbalance, transport, and recombi- 
nation In addition, Butov et al. |1| have observed 

I that, as temperature decreases, the exciton PL fragments 
suddenly into a macroscopically ordered state of nearly 
necklace-like structure: they observed that minimum and 
maximum positions of PL intensity vary nearly period- 
ically along the ordered ring, although the intensity of 
minimum and maximum values are non-periodic. The 
distance between peak positions deviate from the mean 

I period roughly about 20 percent. (For a typical large 
ring observed in Ref. , the distances between peak po- 

I sitions vary roughly between 14 — 22/im.) These random 
variations are not understood yet. 

In this Letter, we derive a nonlinear Gross-Pitaevkii 
equation (GPE) of exciton BEC in the presence of small 
random variation of ring width, and show that the ground 
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state displays a necklace-like structure. Following the ex- 
perimental observation of Butov et al. and Snoke et 
al. lOj we assume that electrons and holes are confined to 
move in a ring. The observed ring radius and width vary 
locally |12|. For simplicity we will model this random- 
ness through a random variation of ring width as excitons 
move along the circumference of the ring with a constant 
radius. This implies that excitons experience a random 
potential as they move along the circumference (smaller 
width gives larger in-plane exciton potential energy). In 
this paper we assume that excitons form a Bose-Einstein 
condensate at temperatures below the quantum degener- 
acy temperature and explore the nature of such a conden- 
sate in the presence of a random potential. We employ 
a mean field approach |^ , which is exact both in the 
low and high density limits. We derive in the low den- 
sity limit a nonlinear Gross-Pitaevskii equation (GPE) 
for excitons and solve it numerically for the ground state 
using the imaginary time evolution method. We find that 
the order parameter of exciton BEC displays a necklace- 
like structure. This is a consequence of a nontrivial in- 
terplay between a repulsive nonlinear term of GPE and 
random potential. Our result suggests that the formation 
of a necklace-like structure observed recently in the PL 
of quantum wells may be stable even in the BEC phase 
at lower temperatures. 

We derive an effective theory of the BEC of two- 
dimensional excitons in a potential trap, which will lead 
to the GPE. We perform a model calculation in which 
the electron-hole interaction to be —g6{x — x'), g > 0. 
We only consider optically active s-wave excitons and ig- 
nore heavy and light hole coupling in the valence band 
[l3l |. We replace the quartic electron- hole interaction in 
the Hamiltonian involving the electron and hole opera- 
tors, -06 (x) and tph (x) , with those containing the pairing 
averages, g{ijjh{y^)'fpe{'x.)), and g(?Aj(x)V'j, (x)) denoted re- 
spectively by the gap functions, A(x) and A*(x). 
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The system is described by the finite- 
temperature Green's functions defined by 
Ge(xr,x'r') = -(rW^e(x, r)VKx', r')), 



Gh(xT,x'r') 



-(r,^^(x,TX(x',T')), ^^C 



-{TrM^,r)M^\r')), and Ft(xr,x'T') 
— (Tt-?/'Jj(x, r)'(/'e(x', t')), where T,- means the time 
ordering in the imaginary time r. These Green's 
functions satisfy 



A(x)* 



V, 



"dr 



1 

2mi 



A(x) 



Ge(xr,x'T') F(xt,x't') 
Ft(xT,xV) -G,,(x't',xt) 



+ (x 
M(x-x')(5(t-t')1 



(1) 



where the gap hmction A(x) = gF{xT^ ,x.t) is obtained 
from the Green's function seh-consistently. Here Ve^hi^) 
is the external potential for the electrons and holes set 
by the potential trap as well as by the random disorder. 
In this model, we also consider a simple situation where 
the effect of the electron-electron and the hole-hole inter- 
actions can be included in Ve,/t(x). Note the presence of 
the magnetic field represented by the vector potential A. 

The coupled differential equations can be written as 
integral equations using the noninteracting Green's func- 
tions, Gg ^ that are obtained for the case where A = 0. 
Iterating the integral equations and keeping up to, say, 
the third order in A, we find 

5-iA*(x)« I d2yQ(x,y)A*(y) (2) 

3 

+ J n'^'y'^^(^'yi'y2,y3)A*(yi)A(y2)A*(y3), 



where Q(x, y) 
and 



iPf^^r' E„ G^(y, X, -c.„)G:!(y, X, 



i?(x, yi, y2, ya) = -mY' E ^^(yi' -^-) 

n 

xG°(yi,y2,w„)G?j(y3,y2,-w„)G°(y3,x,u;„) 

with the fermionic Matsubara frequency uJn- The higher 
order terms in the expansion are given in a similar way. 
The magnitude of the higher order terms will be es- 
timated below. It turns out that it is sufficient to 
keep only up to O(A^) in the present problem. Un- 
like in the derivation of the Ginzburg-Landau equa- 
tions, our GPE is valid in the so-called strong coupling 
regime 16, 17, 18], which is characterized by 

the negative chemical potential /i = /ig + < near 
T = 0, or f3fj, — > — oo. We assume that the effect of 
the magnetic field can be included as a phase factor as 
'^%{^,y:^n) = exp(-i0e,h(x,y))G°^(x,y,w„), where 
'/>e,?.(x,y) = exp(±2|j[A(x) -I- A(y)] • (x - y)) with the 
plus (minus) sign corresponding to the electrons (holes) 
and Gg ^ are the noninteracting Green's function in the 
absence of the magnetic field. This approximation is 



exactly the eikonal approximation made by Gorkov in 
the weak coupling limit 19J . This is reasonable in this 
case also, since A can be considered constant over the 
the size of the exciton. It neglects, however, the terms 
of order A^ [l^, thus not suitable for the strong field 
regime. Within this approximation, the effect of the 
magnetic field disappears in the expressions of Q and 
R in the resulting GPE. Note that, for Ve,h{^) — 0, 
G° f^{x,y,uJn) ~ exp{y/2me,h{lJ'e,h + «w„)|x - y|), which 

decays within the length scale \/It? / {2me,h\fJ-e.h\)- The 
presence of an external potential such a s h. (x) in the 
strong coupling limit was studied in Ref . [l3| ■ As long as 
the external potentials vary slowly over the Bohr radius 
and they are much smaller than the binding energy eo, 
one can replace ^e,h by fie,h ~ (T4,/i(x) + Ve,/i(y))/2 in 
the Green's functions. 

Now, since Q and R are rapidly decaying functions of 
their arguments, we may approximate the first term on 
the right hand side of Eq. I^l as 



A*(x) J d2zg(x,z) + J[v2A*(x)] J rf2zQ(x,z)(x-z)2, 



and the second term as 



A*(x)|A(x)|2 / ndV.i?(x,yi,y2,y3) 



The above integrals are evaluated in the strong coupling 
limit at T = with the assumption fx^^h < 0. In our 
model, excitons move in a Mexican hat potential in addi- 
tion to a random potential represented by Voir). The re- 
sulting GPE is conveniently expressed in terms of the pair 

wavefunction $(x) — y^^A(r) where D = is the 

two-dimensional density of state per spin for the reduced 
mass TO* given by 1/m* = l/rUe + 1/™^- The evaluation 
of higher order coefficients in the GPE is straightforward. 
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Here we present the GPE up to 0($^) as follows: 

+VD{rMr) + ^mr)\H{r) (3) 
— l^frll'^^fr) H \^(r)f^(r) 



where R is the radius of the ring, M — m^+rrih and a de- 
termines the strength of the confinement potential. Note 
that the present GPE is obtained up to 0($'') in contrast 
to the conventional one 0, which contains only 0{^^) 
nonlinear term. The 0(<I>^) term has the opposite sign 
to the 0{^^) term and can become important when ji is 
small. This is, however, not the case in the present prob- 
lem as can be seen below. In general, the magnitude of a 
higher order nonlinear term may be estimated by multi- 
plying the previous lower order term by — ■ The rela- 
tive importance of each nonlinear term will be discussed 
below. The condensate wavefunction is normalized to the 
electron density, given by n(x) = |$(x)p. From the ki- 
netic term of the GPE, we note that the effective mass of 
an exciton is M — + mh, as expected. Note that the 
coefficient of the cubic term in $ is independent of the 
binding energy, which is special to two dimensions. The 
GPE also shows that, in the weak magnetic field regime, 
the exciton does not couple to the vector potential. Mi- 
croscopic parameters appear in the energy /i^, which is 

evaluated to be ^iB = ImKIh^ ' Jd) ^ 1^1 ^"^(r)' 

i_ 

where = e^e 9° is the binding energy of the electron- 
hole pair and eo = with the ultraviolet cutoff A 
used to regularize the delta-function interaction of the 
electrons and holes. This expression for the binding en- 
ergy can be derived from the Schrodinger equation 
for the particle with reduced mass m* and energy —eg 
interacting with the delta function potential of strength 
—g. The strong coupling limit is defined by the chemical 
potential fi = + ^J■h — eo- In the limit fj, — » — eo, 
fiB — eo + /i. This is in contrast with the weak coupling 
BCS limit where /i is approximated by the Fermi energy. 

We rescale the wavefunction $ = ^^E* and the position 
coordinate x = r/a using a length scale a. We define HQ 
by the energy scale associated with the length scale a: 
<j? = The energy scale of the confinement poten- 

tial is fvjj = aa^ . Note that this confinement potential 
alone does not fix the width of the ring since the width 
depends also on the third order coefficient in $ (see be- 
low). Since the number of exciton is / |<I>(r)p(i2r = N , 
the normalization constant is ^ = N^/"^ /a. The number 
of excitons may be estimated from N = {2TrRd)eFDe, 
where c? ~ 1 ^ 10/im is the width of the ring, and 
and ep are the electron density of states and the Fermi 
energy respectively. It is useful to write the GPE in a 
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FIG. 1: Upper panels: Left figure displays the variation of 
|^'(x)p along the circumference of the ring. The dashed line 
shows the random potential yD(x)/50 measured in units of 
hfl. The values of Vd are shifted by an arbitrary amount 
for clarity. The third-order nonlinear coefficient is rii — 10^. 
Right figure shows a two-dimensional contour plot of |^'(x)p. 
Darker regions represent larger values of |'I'(x)|^. The two- 
dimensional position vector is 5 = {xi,X2) measured in units 
of a = 10/im. Lower panels: Same as in the upper panel but 
with rji = 0. Note that the same random potential is used as 
in the upper panels. Only the y-axis scale is difi'erent. As 771 
increases fluctuations of peak values of |^'(x)p are suppressed, 
but the width of the ring increases. 

dimensionless form: 

-iv2vl/(x) + ^rjoix' - xDHix) + yD(x)vE'(x) 

+?7i|^'(x)p^'(x) - 772|«'(x)|%(x) + 773|*(x)|^*(x) 
= /3*(x). (4) 

Here rjo = uj/fl, xq = R/a, and (3 = /is/Sfi. The 
coefficients of the nonlinear terms are ?7i = = 

hft ~ 10~^meV. For exciton densities between 10^ ~ 
10^''cm~2 we estimate rji ~ 10'* ~ 10^. These large val- 
ues will have important consequenses (see below). Since 
\fi\ is comparable to the indirect exciton binding energy 
of few meV we get ep/lfJ^l < 0.01. From these consid- 
erations, we conclude that the nonlinear terms, except 
the third order term, may be neglected. Note that the 
parameter /3, i. e., the binding energy eo is determined 
from the solution of the GPE. 

We now solve the dimensionless GPE for the ground 
state using the imaginary time evolution method. The 
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two-dimensional space is discretized using 512 x 512 grid 
points on a square of area 10a x 10a. The random poten- 
tial is simulated by 250 Gaussian potentials placed ran- 
domly on the ring: Voir) = X^i i where 
Ri are random position vectors on the ring with radius 
R and w — ~ Ifim. The strength of the ran- 
dom potential is taken to a fraction of the strength of 
the confinement potential V — Q.25huj. The correla- 
tion function of the random potential is Vd{x)Vd{x') oc 
_J_^Q-\x~x'\ /2w ^ -y^g gj^^ numerically that the average 
distance between the peaks in |5'(x)p is roughly the cor- 
relation length of the random potential w. The upper 
left panel of Fig. 1 displays how |^'(x)p varies on the cir- 
cumference of ring for rji = 10^ and fku = 2 x lO^hfl. 
A two-dimensional contour plot of |\E'(x)p displays a 
necklace-like structure as shown in the upper right panel 
of Fig. 1 . Positions of modulation maxima are nearly pe- 
riodic while the values of modulation maxima are non- 
periodic. The angular distances between peak positions 
are given roughly as 12 — 20°, which shows a deviation 
from the mean period roughly about 20 percent. These 
results are in qualitative agreement with the experimen- 
tal data (see Fig. 3c of Butov et al §)■ The width of 
the ring increases with increasing coefficient of the non- 
linear repulsive term 771 . In the lower panels of Fig. 1 , the 
solution to GPE is displayed when the repulsive third or- 
der term is absent while the disorder potential is present. 
We observe that the modulation maxima fluctuate much 
more than those of = 10^, and that periodic modulation 
is absent in a two-dimensional contour plot of |5'(x)p. 
Comparing |^'(x)p with the shape of the random poten- 



tial in the left panels of Fig. 1, we can conclude that, as 
the repulsive third order term increases, the peak values 
of |5'(x)p in the local minima of the potential fluctuate 
less from each other. We have also solved the GPE when 
the repulsive third order term is present while the disor- 
der potential is absent. In this case, |^'(x)p is uniform 
along the center of the ring, which implies that the uni- 
form ring state represents a stable BEG of excitons in the 
absence of disorder. 

We have derived a nonlinear Gross-Pitaevkii equation 
(GPE) of exciton BEG in the presence of small ran- 
dom variation of ring width, and have shown that the 
ground state displays a necklace-like structure. Our re- 
sults demonstrate that the interplay between random po- 
tential and the nonlinear repulsive term of the GPE plays 
an important role in understanding necklace-like struc- 
tures. Since the observed fragmented structure appears 
suddenly at low temperature, a collective mechanism of 
quantum origin may be involved, such as the BEG. Since 
we assumed a BEG phase of excitons in our derivation 
our results may be especially relevant in the low temper- 
ature regime of the phase diagram of the ordered state 
proposed by Butov et al . We believe more experiments 
at lower temperature are needed that would elucidate the 
interplay between the non-linear repulsive term and ran- 
dom potential in the formation of the fragmentation. 
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